imshow(p.I);
for i=1:106
    if(p.t(i)==0)
        i;
          hold on;
    plot(p.y(i),p.x(i),'*')
    end
  
end